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ABSTRACT 



We report a glitch-like pulse frequency deviation from the simple spin-down 
^ I law in an anomalous X-ray pulsar (AXP) 4U 0142-1-61 detected by ASCA ob- 

c3 I servations. We also found a significant pulse profile change after the putative 

glitch. The glitch parameters resemble those found in another AXP IRXS 
J170849. 0—400910, in the Vela pulsar, and in other radio pulsars. This sug- 
gests that the radio pulsars and AXPs have the same internal structure and 
glitch mechanism. It must be noted, however, that the pulse frequency anomaly 
can also be explained by a gradual change of the spin-down rate (P) without 
invoking a glitch. 

Subject headings: pulsars: individual (4U 0142-1-61) — stars: neutron — X-ray: 

stars 



^Department of Physics, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8551, 
Japan 



- 2 - 



1. Introduction 

Anomalous X-ray pulsars (AXPs) are a small group of X-ray emitting pulsars (see 
Mereghetti et al. (2002) for a review). They are not likely to be accretion-powered pulsars 
with a binary companion, since the fluxes of the optical counterparts are too small for 
high-mass binary companions or the Doppler modulation due to the binary motion has 
not yet been found. They are also not rotation-powered pulsars, since the observed X-ray 
luminosities (Lx ~ 10^^ — 10^^ ergs s~^) exceed the spin-down energy loss rates of neutron 
stars {E = An'^lP/P^ ~ lO^^.e gj.gg g-i) 

Accretion-powered pulsars without a binary companion (fossil accretion disk) (Chat- 
terjee, Hernquist & Narayan 2000; Alpar 2001) and those with a small binary companion 
(Mereghetti & Stella 1995) have been discussed. However, the pulse timing that originates 
in the accretion disk (large timing noises and/or persistent spin-up periods) is not found. 
AXP IE 2259+586 (Kaspi et al. 2003) and IRXS J170849.0-4000910 (Kaspi, Lackey & 
Chakrabarty 2000; Kaspi & Gavriil 2003) exhibited glitches similar to those of rotation- 
powered pulsars. The timing behavior of AXPs is similar to a solitary neutron star with 
no accretion disk. AXPs are also characterized by soft X-ray spectra (F ^ 2), no radio 
detections, slow rotation periods (5 — 12 s), and burst activities. 

Thompson & Duncan (1996, 2001) proposed the novel hypothesis that SGRs (soft 
gamma repeaters) and AXPs are solitary neutron stars with ultra-strong magnetic fields 
(10^^ — 10^^ G), referred to as "magnetars." X-ray photons are produced in this model by 
the release of the strong magnetic field energy stored in the neutron star crust, or a twisted 
internal magnetic field of neutron stars. Recent discoveries of burst activities (Gavriil, Kaspi 
& Woods 2002; Kaspi et al. 2003), an optical pulsation (Kern & Martin 2002), and cyclotron 
spectral features (Gavriil, Kaspi & Woods 2002; Ibrahim et al. 2002; Ibrahim, Swank & Parke 
2003; Rea et al. 2003) support the magnetar model. Interestingly, the glitch and the pulse 
profile change were observed simultaneously during the outburst of IE 2259-1-586 (Kaspi et 
al. 2003). 

The X-ray source 4U 0142+61 was discovered by Uhuru; its pulsation period of 8.7 s 
was discovered by EXOSAT (Israel Mereghetti & Stella 1994). 4U 0142+61 was one of 
the prototypes of anomalous X-ray pulsars (Hellier 1994; Mereghetti & Stella 1995; van 
Paradijs, Taam & van den Heuvel 1995). A discovery of the optical pulsation supports the 
magnetar scenario by suggesting that the optical emission comes from the magnetosphere of 
the neutron star, not from the accretion disk (Kern & Martin 2002). 

The long-term stability of its intensity, spectrum, and pulse profile were reported through 
a comparison between ASCA observations in 1994 and in 1998 (Paul et al. 2000). The high 
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stability over 4.4 yrs was reported following the long-term monitoring observations of RXTE 
from Nov. 1996 to Apr. 2001 (Gavriil & Kaspi 2002), and the pulse frequency and its 
derivative were precisely determined. However, this RXTE observation was not complete 
due to a 2.0-yr gap from May 1998 to Mar. 2000, dividing its total 4.4-yr span into two 
spans, the "First" (1.3 yr) and the "Second" (1.1 yr) (Tabic 2 in Gavriil & Kaspi (2002) or 
Fig. 1). Gavriil & Kaspi (2002) combined these spans and provided the unified ephemerides, 
"A" and "B," which are qualitatively similar (Table 2 in Gavriil & Kaspi (2002)). 

The ASCA observations in 1998 and 1999 were performed just during the gap between 
the RXTE observations. Those in 1999 in particular were carried out over a period of a 
month, and therefore precise determination of the pulse frequency was possible. We discov- 
ered in our study that the pulse frequency of 4U 0142-1-61 in the 1999 ASCA observations 
was significantly higher than those predicted by the unified ephemerides (A, B) and the first 
ephemeris, and it is also marginally higher than the prediction by the second ephemeris. 
These deviations exclude the unified ephemerides and suggest the presence of a glitch within 
the 2.0-yr gap before the 1999 ASCA observations. In addition, we detected a morphological 
change of the X-ray pulse, in which the pulse profile after the glitch differed significantly 
from those before the glitch. 



2. Observations 

4U 0142-1-61 was observed with ASCA eight times during its mission life. The observa- 
tions were undertaken in 1994, 1998, and 1999. This source was monitored over a period of 
about one month with six exposures during the last term. These observations are summa- 
rized in Table 1. All of the observations were taken by GIS and SIS detectors. Almost all 
GIS observations were taken in the PH mode with 1024 PHA bins and timing bit = at 
high or medium bit rate, in which the time resolutions were 62.5 or 500 ms. We analyzed 
only the GIS-2 and GIS-3 data with high and medium bit rates. 

We obtained all the ASCA observation data for 4U 0142+61 from the High Energy 
Astrophysics Science Archive Research Center (HEASARC) archive site.^ We used the 
SCREENED data, which had been subjected to the standard screening procedure (Revision 
2). We selected the source photons from a circular region with a radius of 5.9 arcmin centered 
on the position with the pulsar peak counts, and we used events with an energy of 0.5 to 10 
keV (44—848 in PI space). The source flux was constant during all observations. The event 
arrival times were corrected to the value at the solar system barycenter using timeconv 



^ftp: / /legacy.gsfc. nasa.gov/FTP / asca/data/rev2 / 
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vl.53 (HEASOFT v5.2). We incremented the times by half of their time resolutions to 
mix the high and medium bit rate data and compensate for the difference of the time 
resolutions, since the times recorded were the leading edge of the time duration of the 
resolution (M. Hirayama; private communication). 



3. Timing Analysis 

We proceeded with the following steps to determine the pulse frequencies (u) and their 
derivatives (z>) for the observation terms of 1994, 1998, and 1999, on the assumption that the 
higher order derivatives were zero. We calculated both values at the reference epochs (to), 
which were selected at the middle of the terms, to minimize the effect of the correlated error 
between u and z>. Hereafter, a subscript "0" denotes a parameter evaluated at the reference 
epoch to. 



3.1. Epoch-folding Method 

The rough {uq, Uq) pairs were estimated by the epoch-folding method, in which the best 
(i/q, ^q) pair was found to maximize the quantity S from each folded light curve: 

s - ±ll^f^ (I) 

{l-2/n)Nj + N/v?' ^ ' 

where n = 16 is the number of bins in one period, Nj {j — 1, 2, ■ • ■ , n) is the number of 
events in the j'-th bin, A'^ = Sj=i total number of events in the observation term, 

and < Nj >— N/n is the mean number of events for the bins. 



3.2. Pulse Arrival Time Analysis 

More accurate determinations of the pulse frequency and its uncertainty were achieved 
by counting the pulse numbers and the pulse arrival times. The pulse numbers 0tot at any 
time t can be expressed by the Taylor expansion around the reference epoch toi as in Gavriil 
& Kaspi (2002). We divided each of the observation terms into short segments. The time 
series was split into short segments of about 30 min duration every ~ 96 min for the terms 
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of 1994 and 1998. Each of the six observations for the term of 1999 (Table 1) corresponds 
to a segment. The pulse profiles of the segments were made by folding the time series with 
the best (z/j, z/j) values of the i-th segment (i = 1, 2, • • • , Useg] Useg- the number of segments) 
obtained by the epoch-folding method cited above, in which the reference "phase zero" epoch 
of each pulse profile was taken at the mid-time of each segment. The template pulse profile 
of the term was generated in the same way using all the segments. Both the template and the 
segment pulse profiles were cross-correlated in the Fourier domains; the Fourier components 
equal to or above the fifth order were discarded. Thus, we obtained the phase offset (AA^j) 
of each segment from the template. The uncertainties of the offsets were estimated by Monte 
Carlo simulation, adding the Poisson errors to the pulse phase bins. 

We were then able to calculate {uq, uq) most accurately by imposing a constraint on the 
segment pulses in which the cumulative pulse phase 0tot(O of the i-th segment relative to 
the template must be an integer (Ni) plus the phase offsets (AA^j). 

The best combination of the Ni{i = 1,2, ■■■ , rigeg) and {uq, z>o) pair was found to minimize 
the following value: 



Here, Ni was selected within the error of 0tot(ii) induced from the uncertainties of {vq, 
z>o) roughly determined by the epoch-folding method. The uncertainties of {vq, z>o) can be 
calculated appropriately, because this analysis uses the fitting method. 



4. Results 

4.1. Pulse Frequency History 

The pulse frequencies and their derivatives obtained by ASCA observations arc listed in 
Table 1. The observations in 1994 and 1998 were too short to determine the z>, and therefore 
we calculated only the v by setting z> = for these terms. The spin-down history is provided 
in Fig. 1. The frequencies in the ASCA observations in 1994 and in 1998 were consistent with 
all of the RXTE ephemerides (Gavriil & Kaspi 2002). However, the frequency for that in 
1999 was significantly higher than the unified ephemerides (A or B) and the First ephemeris, 
and marginally higher than the Second ephemeris of RXTE (Fig. 2) . The frequencies at the 
reference epoch of the 1999 ASCA observation expected by the RXTE ephemerides (Gavriil 
& Kaspi 2002) are listed in Table 2. The same table also reveals the probability of the null 
hypothesis that the frequency of the ASCA observation and those of the RXTE are the same. 
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Table 1: Summary of ASCA Observations 



Start Date 


Exposure 


Reference Epoch 




(UT) 


(ks) 


(MJD) 


(Hz) (Hz s-i) 


1994 Sep 18 (21:35:39) 


18.4 


49614.17288 


0.1151016(8) 


1998 Aug 21 (11:23:48) 


18.9 


51046.69875 


0.1150972(6) 


1999 Jul 29 (16:22:15) 


23.2 






1999 Aug 03 (20:03:29) 


23.5 






1999 Aug 09 (18:37:56) 


11.0 






1999 Aug 12 (22:24:03) 


19.4 


51403.13258 


0.115097645(4) -2.4l°j x 10"^^ 


1999 Aug 21 (19:06:08) 


25.1 






1999 Aug 27 (03:01:19) 


18.1 






Reference epochs are given at the mid-observation. Numbers in parentheses are 1 a uncertainties. 
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Fig. 1. — (Upper panel:) The spin-down history of 4U 0142-1-61. The horizontal and vertical 
axes are the time (MJD) and pulse frequency (Hz). The RXTE "First" and "Second" 
ephemerides are indicated with solid lines. The extrapolations of these ephemerides are 
designated by dotted lines. The ASCA ephemeris in 1999 are represented by a short solid 
line indicated by an arrow. The frequency obtained by the 1994 and 1998 ASCA observations 
are indicated in bars with la errors. (Lower panel:) The glitch magnitude (^) when the 
glitch occurred at the time (MJD) is designated as a la error region. 
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Frequency - 0.115097 (Hz) v.s. TIME (MJD) 

1 1 1 1 1 

RXTE Second 




Time (MJD) 

Fig. 2. — Frequency history of 4U 0142+61 around the 1999 ASCA observation (close-up 
of Fig. 1). The horizontal and vertical axes are the time (MJD) and the frequency minus 
0.115097 (Hz). The spin-down history of 4U 0142+61 expected by the i?XrE observations 
are indicated by the thick lines, where the widths of the lines correspond to la uncertainties. 
The unified ephemerides ("A" and "B") are represented by solid lines: the extrapolations 
from the separate ephemerides ("First" and "Second") arc indicated by the striped lines. 
The error region in the ASCA observation in 1999 is also depicted. 



Ephemeris 


V (Hz) 


a 


Probability 


First span only 


0.115097566(10) 


7.0 


9.7 X IQ-^^ 


Second span only 


0.1150976321(14) 


2.8 


2.4 X 10-3 


Combined "A" 


0.1150976052(12) 


8.8 


< 10-13 


Combined "B" 


0.1150975960(13) 


10.7 


< 10-13 


ASCA 


0.115097645(4) 







Numbers in parentheses are la uncertainties. 



Table 2: Expected frequency at the reference epoch of the 1999 ASCA observation 
(51403.13258 MJD) from the RXTE ephemerides. The significance of the deviation from 
the ASCA ephemeris is also indicated. 
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4.2. 



Pulse Profile 



We searched for the pulse profile changes. We first compared the pulse profiles of every 
segment among the observations in i999 with their template pulse profiles. The pulse phases 
of the segment pulse profiles were aligned with the template using the same cross-correlation 
procedure used in the pulse arrival time analysis (§ 3.2). Each segment pulse profile (D) 
was fitted to the template pulse profile (T) to minimize the reduced statistic by adjusting 
parameters a and b as follows. 



Here, n = 16 is the number of phase bins, and Dj(Tj) is the number of events in the j-th bin 
of the segment (the template pulse profiles). The number of degrees of freedom is u — n — 2. 

There were no significant variations {xl < 1-4 for 14 degrees of freedom, corresponding 
to a probability of > 0.14) of the pulse profiles during the observations in 1999 (Fig. 3). 
However, the template pulse profiles in 1994 and in 1998 differed significantly from that in 
1999 (Fig. 4). The reduced values were 3.8 and 5.9, corresponding to the probabilities 
of 1.8 X 10~^ and 9.3 x 10~^^. The pulse profiles in 1994 and in 1998 were consistent for 



We examined the possibility of whether these pulse profile changes were caused by 
various mixing of the high-bit and medium-bit rate mode data. In fact, the ratio of the 
high-bit to medium-bit rate mode exposures (H/M) ranged from 0.2 to 2.9. We compared 
four pulse profiles created from the 1999 ASCA data with the 1999 ASCA template pulse 
profile to determine how H/M variation affects the pulse profile. These four pulse profiles 
were created from data in which H/M = infinity (H = 0), H/M = 2.8, H/M = 1.2, and H/M 
= 0. The reduced values of the pulse profile differences were 1.8, 1.0, 0.6, and 1.4. This 
investigation verifies that the H/M variation does not significantly affect the pulse profile 
change. 

We also searched for the energy dependence of the pulse profile changes. We divided 
the energy range of 0.5 to 10 keV into a high energy band (3.0 to 10.0 keV) and a low 
energy band (0.5 to 3.0 keV). The threshold of 3.0 keV was determined by the spectrum of 
4U 0142+61. The photon count from the power-law component becomes greater above this 
energy than that from the blackbody component (see Fig. 2 in White et al. (1996)). We 
found that the pulse profile change was only significant in the low energy band. However, 
this does not signify that the pulse profile change only occurred in the low energy band, 
because the photon counts in the high energy band were about 1/10 times lower than that 
in the low energy band. 





x; 
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Pulse Profiles on 1999 ASCA 
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Jul 29 
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Fig. 3.— Pulse profiles of 4U 0142+61 in the 1999 ASCA observations. The vertical and 
horizontal axes are shown in units of the count rates (counts s~^) and pulse phases up 
to 2 periods. (The top panel:) The template pulse profile in the term of the 1999 ASCA 

observation. (The second to the bottom panel:) The segment pulse profiles of the 1999 ASCA 
observations. The reduced vahies (see text) are indicated in the right of the panels. There 
are no significant pulse profile changes during the 1999 ASCA observation. 
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Pulse Profile Difference with Template in 1999 
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Fig. 4. — The template pulse profiles in Sep. 1994 and Aug. 1998 are shown in the upper 
and lower panels. The vertical and horizontal axes are shown in units of the count rates 
(counts s~^) and pulse phases up to 2 periods. The template pulse profile in 1999 is also 
shown by dotted pulse profiles, where the normalizations and baseline offsets were adjusted 
to minimize the reduced (Equation 4). The pulse profiles in 1994 and 1998 differed 
significantly from that in 1999 with xl of 3.8 and 5.9. The residuals from adjusted 1999 
template pulse profiles are also shown. 
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5. Discussion 

We discovered a significant frequency deviation from the simple spin-down law in 4U 
0142+61 for the first time. It is not likely that this deviation was caused by the timing noise, 
because the timing solutions were good during the /^XTi? observations ("First" and "Second" 
spans) and the ASCA 1999 observation. The most likely interpretation of this frequency 
jump is a glitch. The increase of the frequency (Aiy/i/) can be constrained, assuming that 
the ghtch happened between MJD 50893.00 (the last time of the First RXTE observation) 
and MJD 51389.18 (the first time of the 1999 ASCA observation) (Figure 1; lower panel). 
The increase in the spin down rate can be estimated to be Az>/z> = (1.4 ± 0.5) x 10~^ if the 
v difference in the RXTE "First" and "Second" ephemerides (Table 2 in Gavriil & Kaspi 
(2002)) was caused by this glitch. Interestingly, these values are similar to those in a glitch 
of another AXP IRXS J170849.0-400910 (Glitch 1 in Kaspi & Gavriil (2003)). These values 
are therefore also similar to those in the Vela radio pulsar and other young radio pulsars 
(Kaspi, Lackey & Chakrabarty (2000) and references therein). 

We also detected a significant pulse profile change in 4U 0142+61 for the first time, 
which is possibly associated with the glitch. This phenomenon is reminiscent of the out-burst 
that occurred in another AXP IE 2259+586 with regard to the following points (Kaspi et 
al. 2003). (1) The amplitudes of the two peaks in the pulse profile were swapped during 
the outburst in IE 2259+586. (2) IE 2259+586 also underwent a sudden spin-up in the 
magnitude of l^v/v — Ax 10~^. It could thus be interpreted that 4U 0142+61 occasionally 
underwent an outburst as well as a pulse profile change and a glitch just before the 1999 
ASCA observation. Nonetheless, it must be noted that the lack of a pulse profile change 
during the 1999 ASCA observation suggests that survival time of the changed pulse profile 
is more than one month and is longer than the considerably short-lived (~ 6 days) pulse 
profile change observed in IE 2259+586. 

The t\vjv and I\vjv of AXP glitches are similar to those of rotation-powered pulsars, 
and therefore the inner structure of the neutron star and the glitch mechanism may be 
basically the same. Both types of neutron stars may consist of a neutron superfiuid inner 
part and an outer crust part, and the glitch occurs due to the angular momentum transfer 
from the inner superfiuid to the crust by a catastrophic vortex motion. However, the pulse 
periods of AXPs (~ 10 s) are beyond the cutoff of the glitch occurrence, at which periods 
the stress of the Magnus force is not strong enough to trigger the ghtch (see Fig. 8 in 
Ruderman, Zhu &: Chen (1998)). This suggests that there is a difference in the mechanism 
that causes these pulsars to reach the critical state of a glitch; for example, the ordinary 
rotation-powered pulsars reach the critical state due to their spin-down, while the magnetar 
candidates do so due to the stress of their magnetic fields. 
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The pulse profile change suggests two possibihties: (1) Large scale reconnection of the 
magnetic field or deformation of the crust accompanied by deformation of the magnetosphere, 
or (2) Plasma ejection from the crust into the magnetosphere (maybe a burst), which changes 
the distribution of the plasma which emits X-ray photons. The association between the pulse 
profile change and the glitch suggests the following scenarios: (i) The glitch associated with 
the crust cracking (Ruderman, Zhu & Chen 1998; Link & Epstein 1996) or platelet movement 
(Ruderman 1991) causes deformation of the crust, and the pulse profile then changes as 
a result of (1) occurring or (ii) Plasma ejection from the crust associates with both the 
catastrophic vortex motion (a glitch) and the pulse profile changes due to the occurrence of 
(1) or (2). 

We note that the observed ephemeris inconsistency can be explained without assuming 
a glitch. It is possible that the spin-down rate was temporarily smaller for a period be- 
tween the 1998 RXTE observation and the 1999 ASCA observations, as observed in AXP IE 
1048.1—5937 (Gavriil & Kaspi 2004). If the period of the smaller spin-down rate begins at 
the end of the 1998 RXTE observation and ends with the start of the 1999 ASCA observa- 
tion, the spin-down rate in that period becomes minimum. It is z> = —2.49(3) x 10^^^ Hz s~^, 
corresponding to the 5.9 ± 1.2% change from the spin-down rate of the 1998 RXTE observa- 
tion. This lower limit of the i> change is consistent with the finding for AXP IE 1048.1—5937 
(Gavriil & Kaspi 2004). Such behavior can be caused by a change of the magnetospheric 
structure, followed by a torque change and possibly followed by the a pulse profile change. 

6. Conclusion 

An analysis of the ASCA observations suggests the existence of a glitch in 4U 0142+61. 
The ghtch is similar to the ghtch observed in another AXP IRXS J170849. 0-400910 and 
other young radio pulsars. This indicates that the radio pulsars and AXPs have the same 
internal structures and glitch mechanisms. However, the source of the stress triggering the 
glitch may be different. We also discovered a pulse profile change, possibly associated with 
this glitch. This suggests that AXPs have an association between the internal structure and 
the magnetosphere that does not exist in radio pulsars. 

This work was partially supported by a 21st Century COE Program at TokyoTech 
"Nanometer-Scale Quantum Physics" by the Ministry of Education, Culture, Sports, Science 
and Technology. This research was also supported in part by Grant-in- Aid for Scientific 
Research(C) (15540239). 
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